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Abstract. We present the numerical implementation of run-time injection of Cosmic 
Rays energy, their spatial advection and their dynamical feedack on baryonic gas in the 
cosmological grid code ENZO. We discuss the results of its application to large-scale 
simulations showing that the CR energy inside clusters of galaxies is small compared 
to the gas energy (less than a few percent), while the ratio is larger near the accretion 
regions of clusters and filaments (~ 0.1 - 0.3). CR feedback has a small, but significant 
impact on the X-ray emission and Sunyaev-Zeldovich eff'ect from clusters. 



1. Introduction 

The detection of large-scale non-thermal emissions in galaxy clusters proves the exis- 
tence of magnetic fields and relativistic particles (CR) in the intra cluster medium (e.g. 
Ferrari et al. Merger events in clusters drive shock waves in the intra cluster 



medium, and this can inject primary CR particles and may oflFer a simple explanation 
of the observe d connection between large-scale radio sources with arc-like shapes (e.g. 



"radio relics", Briiggen e t all 120111 for a recent review) and the underlying disturbed 
cluster morphology. The mechanism to convert kinetic energy from the ~ Mpc scale 
of shock waves into the acceleration of CRs may be ex plained by the DiflFusive Shock 
Acceleration (DSA, e.g. "BeH""l978'; "KangAJQnes|l2QQ7). According to this model, a 
pool of supra-thermal particles is produced by collision-less shocks in rarefied space 
plasmas, and are accelerated to higher energies by the interactions with Alfven waves 
in the post-shock region. We present our implementation of a model of DSA in the 
public release of the cosmic code ENZO 1.5 (see Sect]!]), simulating the impact of CR 
physics in the evolution of large-scale structures (see Sect|3]). 



2. Numerical methods 

All simulations presented in this w ork were performed using the grid-based and adap- 
tive mesh refinement code ENZO (iNorman et al.ll2007l : [Collins et all 1201 d ). ENZO is 
currently developed developed by the Laboratory for Computational Astrophysics at 
the University of California in San Diego (http://lca.ucsd.edu ). In our work, we im- 
plemented the treatment of CR energy injected at shocks by adding several routines, 
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Figure 1 . Left and central panel: 1-D shock tube test for an initial mixture of gas 
and CRs (shown ar e density, gas, CRs and total pressure). The red lines show the 
exact solution from iMiniatil (l2QQ7h . Right panel: average pressure ratio Pcr/Pg as a 
function of gas overdensity for a cosmological run at fixed grid resolution (Ax =156 
kpc/h) for diff'erent cosmic epochs. Additionally shown is dotted-red is the result at 
z=0 for an AMR resimulation of a cluster, with a maximum resolution of Ax = 25 
kpc/h. 



physical fields and new sets of e quations to the public 1.5 release of ENZO, following 
a two-fluid model approach (e.g. iKang & Jonesll2007l ). 

An on-the-fly shock finder was developed to detect shocks in simulations at each 
time step. This shock finder is ba sed on a 3-D analysis of gas pressure, similar to 
other methods in the literature(e.g. IRvu et al.]l2003l : ISkillman et al.ll2008l : IVazza et aT 



20091 ). Once the Mach number, M, is measured, we estimate the total energy flux of 



accelerated CR protons as: 0cr = t](M) • PuVs/2 (pu is the pre-shock gas density, Vs 
is the shock velocity , and the function 7](M) models the eflftciency of injection, as in 
Kang & JonesI (l2007l )). We inject CR energy in each post-shock cell as Ecr = 0crA/Ax^, 



where A^ is the time step and Ax is the cell size; a corresponding amount of energy 
is also decreased in the post-shock region to conserve total energy. The CR energy is 
advected with the flow assuming no diflFusion and no energy losses for CRs, using an 
adiabatic index of ycr = 4/3. The total pressure P = Pg + Per (where Per = (7cr - ^)Ecr 
and Pg is the gas pressure) is fed into the PPM Riemann solver of ENZO (rather than 
the gas pressure). The composite fluid gas-hCR in the simulation obeys an eflFective 
adiabatic index, yeff < y = 5/3. To preserve the stability of the method, we also enforce 
an additional condition on the time step (which avoids too strong a modification of the 
post-shock during a single time step). For a full description of the adopted methods, 
and of their tests against analytical solutions we refer the reader to a forthcoming paper 
(Vazza et al., in prep). In FiglH we show the comparison of a 1-D shock tube test and 
the numerical solution obtained with the Glimm-Godunov's method by lMiniatil (l2007l ). 
showing the very good performance of our new scheme. 



3. Results 



We present here the results of non-radiative cosmological runs for a concordance" 
ACDM cosmology, with parameters Qdm = 0.226, = 0.044 and Qa = 0.73; the 
normalization of the primordial index of density fluctuations was se t to erg = 0.8. In 
our runs, we used an acceleration efl&ciency of iKang & JonesI (l2007l ): the injection of 
CR energy is modelled since the start of reionization in the simulated volume (our tests 



Cosmic Rays in large-scale structures 



3 



however suggest that the initial epoch for the injection of CR is not an important pa- 
rameters, provided that Zcr » !)• In figure[2](top panels), we report maps of the spatial 
distribution of gas and CR pressure for the simulated volume of 80 Mpc at the uniform 
resolution of 156 kpc/h at z=0 (512^ cells). Despite the overall morphological similar- 
ity between thermal gas and CR structures, the CR energy is much more concentrated 
around large-scale structures than gas energy; this is explained because CRs are injected 
only in the downstream region of shocks, and because their adiabatic compression dur- 
ing the infall onto collapsed structures is smaller than that of gas, due to their softer 
equation of state. The average pressure ratio between CRs and gas within the simu- 
lated volume has a maximum of Pcr/Pg ~ 0.2 - 0.3 at the cosmic gas critical density, 
and it declines smoothly towards the centre of massive galaxy clusters, approaching 
Pcr/Pg ~ 0.05 - 0.1 there (see right panel in Figil]). Re- simulations of individual clus- 
ters employing adaptive mesh refinement ( with maximum resolution of ~ 25 kpc/h) 
provide an even smaller pressure ratio for the core of clusters, Pcr/Pg ~ 0.01 - 0.03 
(Fig|2]and right panel of Fig{T])- The dynamical role of CR energy is negligible for most 
of the simulated volume; however, the final distribution of yeff as a function of radius 
from the center of structures is not uniform, and decreases towards the outer regions, 
making the accretion regions of clusters more compressible. This systematically pro- 
duces a slightly more relaxed matter distribution within clusters, compared to standard 
runs without CR injection, and this has a small but significant (~ 5 - 10 per cent) impact 
on the measured X-ray emission profiles and Sunyaev-Zeldovich effect (see Vazza et 
al. submitted). 



4. Conclusions 

We presented a numerical development in the grid AMR code ENZO 1.5, which allows 
us to treat at run-time the injection of CR energy at shocks, their spatial advection and 
their dynamical (pressure) feedback on the evolution of the baryonic gas. This approach 
provides an important and complementary approach to similar techniqu es applied to 
Smoothed Particles Hydrodynamics methods (e.g. Pfr ommer et al.|[200fl ). Simulations 
employing fixed resolution or adaptive mesh refinement have been run, leading to the 
detection of small but non-negligible effects respect to standard large-scale structures 
simulated without CR physics. The application of these techniques to the problem of 
large-scale non-thermal emission (as for example in radio relics) is underway. 
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Figure 2. Maps of gas and CR pressure for the simulated volume with the side 
of 80 Mpc (top panels) and for the AMR resimulation of a galaxy clusters (lower 
panels, side of 10 Mpc). The color coding is in log[arb.units]). 
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